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Abstract 

We consider the mean distribution functions $(r|£), <J>B( r |^)> and §§{r\(), 
giving the probability that two sites on the incipient percolation cluster, on 
its backbone and on its skeleton, respectively, connected by a shortest path of 
length I are separated by an Euclidean distance r. Following a scaling argu- 
ment due to de Gennes for self- avoiding walks, we derive analytical expressions 
for the exponents g\ = df + d m i n — d and gf = gf = 3d m i n — d, which determine 
the scaling behavior of the distribution functions in the limit x = r ji v <C 1, 
i.e. <f>(r\£) cx f m , <f> B (r\£) oc i~ M , and $g(r\£) oc £~ M x^, with 
v = l/d m { n , where df and d m \ n are the fractal dimensions of the percolation 
cluster and the shortest path, respectively. The theoretical predictions for g\, 
gf and gf are in very good agreement with our numerical results. 

PACS numbers: 05.20.-y, 64.60.-i, 05.40.+j 
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Percolation constitutes a useful model for a variety of disordered systems in many fields of 
science displaying both structural disorder and self-similarity (i.e. fractal behavior) within 
some range of length scales In many circumstances, the knowledge of the internal 
structure of percolation clusters is required, as for instance in the study of transport processes 
near the percolation threshold p c , where the complex topology of the available conducting 
paths play a crucial role 0-|5[]. 

It is known that at the percolation threshold p c , the incipient infinite cluster displays 
fractal behavior over all length scales, i.e. its mass s contained within a distance r from a 
given cluster site chosen as the origin, averaged over many origins, scales as s oc r df , where 
df = 91/48 in two dimensions, df = 2.524 ± 0.008 in three dimensions, and df — 4 above 
the critical dimension, i.e. when d > d c = 6 [Q. A second, useful metric is the "chemical" 
distance £ between two cluster sites , defined as the length of the shortest path connecting 
them. It is found that the mean distance r between two cluster sites, averaged over many 
pairs of sites, behaves function of £ as r oc 

£iAUn ; w h ere d min = 1.130 ± 0.004 in d = 2 
II j ^min = 1-374 ± 0.004 in d = 3 0, and d min = 2 when d > d c , is the so-called fractal 
dimension of the shortest path. From the above scaling relations follow that in "chemical" 
space, the mass of the cluster scales with distance £ as s oc £ de , where di = df/d m i n , with 
d e = 2 when d > d c [[J . 

The incipient infinite cluster exhibits a variety of substructures that are self-similar as 
well [0. A prominent example is the backbone of the cluster, defined as the subset of 
cluster sites that can carry a current when a potential difference is applied between two 
sites (see || and references therein). Thus, the structure of the backbone alone determines 
the conductivity of the whole percolation network between two sites. The structural and 
dynamical properties of the backbone of the incipient cluster have been studied recently 0. 
A second cluster substructure, denoted as the skeleton (a subset of the backbone, also called 
the "elastic" backbone) is defined as the union of all shortest paths between the two cluster 
sites. 

In this Rapid Communication, we extend our previous studies of the structural properties 
of the incipient infinite cluster [jHJ and its backbone || in two and three dimensions. We 
consider the structural distribution function $(r|£) for the incipient infinite cluster, where 
$(r|£) dr is the probability that two cluster sites connected by a shortest path of length £ 
are at Euclidean distance between r and r + dr from each other in space. The probability 
distribution $>(r\£) is normalized according to J r d_1 $(r|£) dr = 1, and is found to obey an 
scaling behavior with the variable x = r/£ u of the form <&(r\£) = £~ ud f(x) (see e.g. P,|IT)|,|Tl"|]), 
where v = l/d min . Here, we draw our attention to the limit i < 1, where the scaling function 
f(x) follows a simple power law, f(x) oc x 91 , i.e. 

1 / r \ gi 

*(r\£) ex ^ i^j , for r/r « 1 . (1) 

Similar scaling forms for the substructural distribution functions as a function of x = 
r/t, <&B{r\£) = £~ M /b(^) for the backbone and $s( r l^) = fs{ x ) f° r the skeleton, are 
expected J|||. In the case x C 1, the corresponding scaling functions, /b(^) and fs{x), are 
found to behave as /b(^) oc x 91 and fs(x) oc x 9 *, respectively, yielding 

W)^(^) Sl and $ s (r|£) oc ^ {^f , for r/t « 1 . (2) 
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Numerical results (see Refs. [§,fjlj and below) indicate that gi < gf = gf in both two and 
three dimensions. For d > d c , one expects the mean field (MF) values gi = gf = gf = 0, since 
percolation clusters behave similarly to simple random walks above the critical dimension 

We first study the above defined distribution functions numerically, both in two and 
three dimensions. To this end, we generate large percolation cluster at p c on square and s.c. 
lattices, respectively, using the well-known Leath algorithm |T2|]. To identify the backbone 



and skeleton of the cluster, we apply an improved version || of the well-known burning 
algorithm ||. We perform averages over more than 10 5 clusters, which are grown until they 
reach a maximum of chemical shells £ max = 2000 in d = 2 and £ max = 1000 in d = 3. The 
results for $(r|£), $b(j|^), and $s(V|^) are shown in Figs. |], §, and ||, respectively. For the 
incipient infinite cluster we obtain g\ = 1.04 ± 0.05 in d = 2 and g\ = 0.88 ± 0.05 in d = 3 
(see also 0). For the backbone, we find gf = 1.34 ± 0.10 in d = 2 and gf = 1.08 ± 0.10 
in d = 3 (see also [[J). In addition, our results suggest that $b(?"|^) an d $s( r l^) coincide, 
within the accuracy of the present data, and as a result, the values of gf for the skeleton are 
indistinguishable from those of the backbone, i.e. gf = gf. These results are summarized in 
Table I. 

To estimate values for the exponents gi, gf and gf analytically, we follow a method 
similar to the one discussed by de Gennes |13[] for determining the structure of self-avoiding 
walks (SAW) of iV steps. The latter is described by the probability distribution function 
-PsawM^O = N~ ud /sAw(y)> with y = r/N u , where PsawM-^O dr gives the probability that 
the two end points of a SAW of fixed length N (i.e. the first and the N + 1 monomers) are at 
a distance between r and r + dr. Here, v is the Flory exponent, v = (d + 2)/3 for d < 4 and 
v = z/ MF = | for d > 4, and /saw(?/) is the scaling function, with /saw(z/) ^ V 9 when j/< 1. 
For SAW defined on the lattice, de Gennes argues that the behavior of /sAw(y) for y <C 1 
can be obtained by considering the probability -PsAw( r "~ > that a SAW of iV ^> 1 steps 
returns close to its starting point (origin), which can be written as 

PsMr - m k SS! ■ foriv>>1 ' (3) 

where N^^(N) oc N~ ud z N is the number of SAW of length iV returning close to the origin 
and Nsaw(N) oc iV 7-1 z N is the total number of SAW of length N. Here, z is the effective 
coordination number of the lattice, and 7 is the enhancement exponent, with 7 = 7mf — 1 
for d > 4. 

As noted by de Gennes ||13[| , the enhancement factor iV 7-1 occurs only in the denominator 
of the ratio N^J^(N) / N S aw(N) , but not in the numerator, indicating the "difficulty" for a 
SAW to return near to its starting point. Note that this missing enhancement factor in the 
numerator can be viewed as corresponding to its mean-field value, 7V 7MF-1 = 1, and one can 
write equivalently 

1 7V7MF-1 

Wr-.llJVja^^, foriV»l, (4) 

corresponding to the behavior /sAw(y) oc y 9 , with g = (7 — l)/v. This observation suggested 
us a procedure for describing the structural function of the incipient percolation cluster and 
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its substructures analytically, in the case r/£ u <C 1. We consider the incipient percolation 
cluster first. 

Let us generalize Eq. (B) to percolation clusters by writing the distribution function 
$(r|£), for a chemical distance £ ^> 1 and Euclidean distance r — > 1, as 

*( r ^M«^$T> for£>>1 ' ^ 

where U(£) plays the role of the function A r7_1 in Eq. (||), and IImf(^) denotes its mean 
field value. Here we argue that, to a first approximation, II (£) is given by the probability 
that the two chosen sites are on a cluster of chemical size £. Therefore we relate H(£) to the 
probability distribution of cluster sizes sn(s), which is known to behave as sn(s) oc s _ ( T_1 \ 
with t = 1 + d/df for d < d c , and t mf = 5/2 |L|. Hence, n(£) is given by n(£) oc sn(s) ds/d£, 
and noting that s oc £ dt , we find U(£) oc £~Mr-2)-i for d < ^ and n MF (Q oc £~ 2 . Thus, 
Eq. (|5p becomes 

^- 1 ^«^P5^=r«p5^ (T - a) - 1 . fo ^» 1 - ( 6 ) 

Comparing this result with the one obtained from Eq. (|l|) in the limit r — > 1, yields — P<7i = 
d,(r - 2) - 1, i.e. 

gi = d f + d min - d , (7) 

which predicts #i = 1.026 ± 0.004 for d = 2 and #i = 0.898 ± 0.008 for d = 3. These 
theoretical values for g\ are in very good agreement with our numerical results (cf. Fig. [I] 
and Table |). Note that Eq. (|7|) yields by construction g x = for d > d c , as required. 

The above argument can be applied straightforwardly to the backbone and the skeleton 
of the incipient cluster, where now analogous equations to Eq. (|5|) can be written for $B( r ~~ > 
1\£) and $s( r — * with H(£) replaced by n B (^) and n s (£), respectively. In the case of 
the backbone, we argue that n B (^) oc n(s) ds/d£, with n(s) oc s~ T , and s oc £ dl as for the 
incipient cluster. Note the absence of the factor s in the expression for Hb(^), reflecting 
the fact that the backbone represents a subset of the incipient cluster having a vanishing 
measure when s — > oo [|TJ[]. Since the same argument applies to the skeleton, we have that 
U s (£) = Il B (£), yielding 

$ s ( r -> 1|£) S $ B ( r -> 1|£) , for£>l, (8) 

in agreement with the numerical results shown in Figs. [2] and |3|. In terms of Hb(^), $b(^K) 
in the limit r — > 1 is given by 

* B(r - m K wWr ■ for<>>1 - (9) 



and with Hb^) oc 



(r-l)-l 



$ B(r -> ex ^ ^"-D-i « ^ (T ' 1} ' 3 ' for£>>1 - (10) 
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Comparing this result with the scaling form for $b(V|£) given in Eq. (|2|) in the limit r — > 1, 
yields — vgf = dfij — 1) — 3, i.e. 

gf = 3d min - d , (11) 



B 
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predicting gf = 1.390 ± 0.004 in d = 2 and gf = 1.122 ± 0.004 in d = 3, with gf = g[ 
remarkable agreement with the numerical results (cf. Figs. [2] and [3], and Table [I]). Note also 
that from Eqs. ([□]) and (H) one obtains gf = gf = for d > d c , as expected. 

In summary, we derive the analytical expressions g\ = df + d m i n — d and gf = gf = 3rf min — 
d describing the scaling behavior of the structural distribution functions, $(r|£) oc £~ ud x 91 , 
$b(^|^) oc i~ vd x 91 , and $s( r l^) i~ vd x 91 , of the incipient percolation cluster, its backbone 
and skeleton, respectively, at the critical concentration p c in the limit i; e r/f <C 1. Here, 
v = l/rf min , and df and d min are the fractal dimensions of the incipient percolation cluster 
and the shortest path, respectively. Note that from the above expressions for the exponents 
gi, gf, and gf follow that the corresponding distribution functions for £ 3> 1, in the limit 
r -> 1, scale as $(r -> 1\£) oc £~^ +1 ) and $ B (r = $ s (r 1\£) oc £~ 3 , the latter 

being independent of the lattice dimension d. We note that the result $(r — > 1\£) oc 
for £ ^> 1, based on numerical simulations, was also suggested for two other variants of 



percolation, invasion percolation with as well as without trapping and seems therefore 
to be more general. 
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TABLES 







d = 2 




d = 3 


d = 6 




simulation 


theory 


simulation 


theory 


exact 


91 


1.04 ±0.05 


1.026 ±0.004 


0.88 ±0.05 


0.898 ± 0.008 





9f 


1.34 ±0.10 


1.390 ±0.012 


1.08 ±0.10 


1.122 ±0.012 





9t 


1.34 ±0.10 


1.390 ±0.012 


1.08 ±0.10 


1.122 ±0.012 






TABLE I. Summary of the values for the exponents g±, gf and gf obtained from the numerical 



simulations and the analytic expressions derived in the text. 
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FIGURES 
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FIG. 1. Scaling plot of the probability distribution function i ud 5>(r|£) vs r j£ v for the incipient 
infinite cluster in the following cases: (a) d = 2, i = 1000 (circle), £ = 1400 (diamond), and 
I = 1800 (square), and, (b) d = 3, £ = 400 (circle), £ = 600 (diamond), and £ = 800 (square). The 
plots are based on averages over more than 10 5 cluster configurations, for clusters grown up to a 
maximum chemical distance £ max = 2000 on a square lattice (d = 2) and £ max = 1000 on a s.c. 
lattice (d = 3). The straight lines represent our fits for £ ud <E>(r|£) = f(x) when x = r jl v <^ 1, and 
have the slopes g\ = 1.04 in (a), and gi = 0.88 in (b). 
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FIG. 3. Same as in Fig. [I] for the probability distribution function °f the skeleton of 

the incipient cluster. The straight lines have the slopes gf = 1.34 in (a), and gf = 1.08 in (b). 
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